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Abstract 

We consider the expectation value of a local operator on a strip with non-trivial bound- 
aries in 1+1 dimensional massive integrable QFT. Using finite volume regularisation in 
the crossed channel and extending the boundary state formalism to the finite volume 
case we give a series expansion for the one-point function in terms of the exact form 
factors of the theory. The truncated series is compared with the numerical results of 
the truncated conformal space approach in the scaling Lee^ang model. We discuss the 
relevance of our results to quantum quench problems. 
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1 Introduction 

Finite size effects play a central role in quantum field theory and statistical physics. Apart from hav- 
ing direct relevance to statistical physics models in finite volume, describing for example boundary 
critical phenomena or percolation problems, they naturally appear in the description of systems at 
finite temperature. For example, two-dimensional Euclidean field theories with a finite, periodic di- 
rection provide a framework for studying one-dimensional theories at finite temperature. Moreover, 
boundary phenomena and finite volume systems can play an important role in the understanding of 
quantum quenches: in certain cases the boundaries play the role of the initial and final states of the 
non-equilibrium problem [1, 2, 3]. In addition, in many cases, especially for numerical simulations, 
the system under consideration is put in a finite volume box. In this case it is essential to understand 
the finite size behaviour of various quantities. 
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Correlation functions are very important both at finite and infinite volume because they encode 
a lot of non-trivial information about the spectrum and the interactions in the theory. Calculating 
correlators is not a simple task even in integrable theories. In 1+1 dimensional integrable quantum 
field theories the form factor approach provides an efficient way to calculate correlation functions 
at zero temperature and in infinite volume. The essence of the method is inserting a complete set 
of asymptotic states into the correlators and then making use of the explicit forms of the appearing 
matrix elements of the operators, the so-called form factors, which are known in many integrable 
models. 

However, for finite temperature or equivalently, for finite geometry the applicability of the form 
factor approach at present is somewhat constrained. There has been considerable progress in free 
theories (for example the quantum Ising model), where the non- interacting nature of the theory 
allows for calculating correlation functions both in the finite temperature setting [4] or in the 
presence of a boundary [5, 6]. However, interacting theories pose technical and conceptual difficulties, 
which have not yet been overcome. 

In the finite temperature case one particular approach was developed by LeClair and Mussardo 
[7]. They proposed an integral series for correlation functions based on the exact form factors and 
the Thermodynamic Bethe Ansatz equations. Their result for one-point functions were checked in 
particular examples [8, 9, 10], and then a highly non-trivial check was given by confirming it up 
to the third order in the low-temperature expansion using finite volume regularisation [11]. For 
two-point functions, however, some counterexamples were found where the formalism does not seem 
to work [8, 9] and the problem is still far from being settled. 

In this work we address the generalised problem of one-point functions in finite volume where 
instead of periodic boundary conditions we consider non-trivial boundary conditions, thus we have 
a strip geometry instead of a cylinder. In this sense our work can be a first step in finding an 
expression similar to the LeClair-Mussardo formula in the boundary case. Expectation values in 
the presence of boundaries in general depend on the position of the operator in question^. The exact 
determination of them is thus non-trivial; for one boundary the problem has comparable difficulty 
to that of the two-point functions in infinite volume; and similarly, in the presence of two boundaries 
the technical problems resemble the case of a three-point function or a two-point function at finite 
temperature. 

In the work [13] the authors considered the case of one boundary, or in other words the prob- 
lem when the operator is much closer to one of the boundaries than the volume of the system. 
They developed a form factor expansion which was later used in [14, 2] to study certain prob- 
lems in condensed-matter systems. In addition, the analytic results were confronted in [13] with 
the numerical Truncated Conformal Space Approach (TCSA). Our expansion reproduces and goes 
beyond their result, and we will show that at small enough volume our results give a remarkable 
improvement. Our numerical comparison with the TCSA results is very similar to their method. On 
the other hand, the theoretical determination of the spectral series is more involved: there appear 
conceptual (and also technical) difficulties, which are not present in [13]. Most importantly, one has 
to deal with certain singularities of the form factors and also the divergent contributions to the 
partition function. In this work we make use of a finite volume regularisation scheme, which first 
was used in interacting field theory in [15, 11] and later in [16, 17]. 

As remarked above, integrable boundary field theory can be used to investigate quench problems. 
Recent papers [2, 3] consider the time evolution of one-point functions after certain types of global 
quenches. As we will show, their problem and method have similarities to ours and thus our results 
also have relevance to quench problems. 

^This has to be contrasted with the problem of expectation values of boundary operators considered for example 
in [12]. 
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The outline of the paper is the foUowing. In section 2 we introduce the finite volume boundary 
states and as a first application we give a series expansion of the finite volume ground state energy. 
In section 3 we turn to the problem of one-point functions. After describing the method in detail 
we present the explicit calculations of the first terms in the form factor expansion. At the end of 
the section we discuss the relation of our expression with existing results in the literature and we 
also explain its connection to quench problems. Section 4 is devoted to the numerical comparison 
between the spectral expansion and the TCSA data for the scaling Lee-Yang model, considered 
as a perturbed conformal field theory. First we collect the relevant properties and formulae of the 
model, then we discuss briefiy the TCSA method and how one-point functions can be obtained in 
this framework. Finally we compare the form factor and TCSA results for various combinations of 
boundary conditions and strip widths. Our conclusions are given in section 5. For some technical 
details and for a collection of our final result the reader is referred to the appendices. 



2 Boundary states and expectation values 

For the sake of simplicity let us consider an integrable relativistic field theory with only one particle 
species with mass m and two-particle S-matrix S{9). We are interested in the case where the theory 
lives on a finite line segment of length R with boundary conditions at the edges that do not spoil 
the integrability. 

If B denotes such an integrable boundary then the scattering of an incoming multi-particle 
state'^, 

\Ai9i)A{92)...A{en))B , (2.1) 

is purely elastic, i.e. the set of outgoing rapidities is {— ^i, —62, ■ ■ ■ , —9n} and the scattering can be 
described by a product of one-particle refiection amplitudes which are defined as [18] 

\A{d))B = RB{0)\A{-e))B . (2.2) 

The main objective of this paper is to evaluate the vacuum expectation value 

{Oixr/ (2.3) 

in the presence of the integrable boundaries a and /3 where x € [0, R]. Here the expectation value is 
taken with respect to the vacuum state of the finite volume Hamiltonian H"^^ (see fig. 1(a)). In a 
more general setting one can also take excited states corresponding to boundary bound states but 
we do not elaborate on this case here. 

One can consider the same quantity (2.3) after a Euclidean rotation. In this picture R plays the 
role of the Euclidean time variable and the expectation value is given by the formal expression 

^^^'''^R - {B^\e-^R\Bp) ' ^^■^> 

where H is the infinite volume Hamiltonian and \Ba) and \Bp) are the boundary states correspond- 
ing to the boundary conditions a and /?, respectively. When they do not contain zero-momentum 
particles, they can be expanded in the asymptotic multi-particle basis as [18] 

= AA.exp (I ^K,{e)A{-e)A{e)^ |0) , (2.5) 



'The subscript B indicates that the state is in the Hilbert space of the boundary theory. 
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Figure 1: One-point function in the original strip geometry (a) and the finite volume regularisation 
performed in the crossed channel (b). The upward arrows denote the (imaginary) time evolution 
generated by the corresponding Hamiltonians. 



where 

K,{e) = R,{i7r/2-9) , j = a,P 
and A{9) are the Faddeev-Zamolodchikov creation operators satisfying the commutation relations 

A{9i)A{e2) = s{ei - e2)A{e2)A{e,) . 

The integrals run from — oo to oo unless otherwise stated. The amplitudes Ki[9) satisfy the "bound- 
ary cross-unitary equation" [18] 

K^{e) = s{2e)Kj{-e), 

which serves as a consistency relation of (2.5). The behaviour under complex conjugation is given 

hy Kj{er = K,{-e). 

In general there may be additional contributions to the boundary state involving zero-momentum 
particles, which can be associated with a non-zero coupling of a single particle state to the boundary 
in the original channel. In this case there is a corresponding pole in the reflection factor at = m/2: 



Rj{e)r^- 



29-^/2 



and the boundary state is given by 



|i?,)=AG exp (5,-4(0) + j ^K, 



ie)Ai-e)Aie) \o) , 



(2.6) 



where gj represents the one-particle coupling to the boundary. Originally it was argued that gj = gj 
[18], however it was found numerically in [13] that the proper normalisation is gj = gj/2. This claim 
was later proven on general grounds in [19, 20]. The normalisation constants Mj are infinite in the 
infinite volume system. However, this is not a serious problem as they drop out from the calculation 
of the vacuum expectation values and one may set them to unity. 

The general strategy to evaluate (2.4) now is the following. We substitute expression (2.6) for 
the boundary states and expand the exponentials. The resulting multi-particle states are eigenstates 
of the Hamiltonian, and the matrix elements of the local operator O between these states are the 
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form factors, which can be determined in principle and which are actually known in many integrable 
models. 

The special case for the i? — )■ oo limit of (2.4) with x fixed (one-point function in the presence 
of only one boundary) was investigated in [13]. In this case there are no conceptual and technical 
difficulties and the spectral expansion described above (valid for = 0) reads 

(Oi-))=fl^^ I ^...|^i^2''n(-^l,ei,...,-^n,^n)n(^^«(^0e-^'"^°^''") > (2-7) 
n=0 ■ i=l 

where 

F^{ei,...,em) = (o|o|0i,...,M 

are the elementary form factors of O in infinite volume. In integrable models they are meromorphic 
functions of the rapidities which satisfy the so-called form factor axioms (Watson equations), which 
can be considered as axioms for the form factor bootstrap. Supplied with the principles of maximum 
analyticity and the cluster property they contain enough information to determine the form factors 
completely. In a theory with only one particle species the form factor axioms read: 

I. Lorentz invariance: 

F^iOi +A,...,9n + A) = e'^F^ie,, ...,0n), (2.8) 

where s is the Lorentz spin of the operator O. In this work we only consider scalar operators 
corresponding to s = 0. 

II. Exchange: 

F^{6i, . . . , Ok, Ok+i, . . . ,9n) = S{0k - dk+ijF^ {6i, . . . , 6k+i, Ok, ■ ■ ■ , 6n) ■ (2.9) 

III. Cyclic permutation: 

FO{9i + 2i7r,e2,...,en) = F^{92,...,en,9i) . (2.10) 

IV. Kinematical singularity 

- i Res F^ {0 + i7r,e' ,01,. . . ,en) = (l -f] 3(6 - ek)] F^{e^,. . . , On) . (2.11) 

V. Dynamical singularity 

- iRes F^^^{e + iu,e' - iu,ei, . . . ,en) = TF^^^{e,ei, . . . ,en) , (2.12) 

e=e' 

where F is the on-shell three-particle coupling corresponding to a bound state pole of the S-matrix 

Sie ~ i2u) , 6* ~ 2m . (2.13) 

— i2u 

In a theory with only one particle, like the scaling Lee^ang model the only possibility \s u = 7r/3. 

We also note that all form factors can be expressed in terms of the elementary form factors with 
the help of the crossing relation 

-^mn(^l' • • • ' ^ml^l' • • • ' ^n) = -^m-l,n+l (^1 ' • • • ' ^m-ll^m + "IT^ , 6l, ■ ■ ■ , On) 

n k—1 

+ (27r(5(6'^ - Ok) S{Oi - Ok) x F^_i^„_i(6'i, . . . , 0^_i\Oi, . . . ,0k-i,0k+i ■ ■ ■,0n)^ ■ 

k=l 1=1 



6 



For an introduction to the form factor bootstrap program we refer the reader to the review of 
Smirnov [21] and to papers [22, 23, 24, 25]. 

In the simultaneous presence of the two boundaries the expression (2.4) is ih-defined. There 
are divergent contributions both to the numerator and the denominator, which are associated with 
disconnected terms of the form factors and with various contributions to the partition function. One 
way to evaluate (2.4) is to introduce finite volume as a regulator of the singular contributions. This 
method was successfully used in [11] to derive a rigorous low-temperature expansion for one-point 
functions at finite temperature, which in the crossed channel corresponds to vacuum expectation 
values in a finite volume with periodic boundary conditions. In this sense the present work is an 
extension of the approach of [11] to the more general setting of arbitrary (integrable) boundary 
conditions. 

We introduce a large finite volume L in the crossed channel (with periodic boundary conditions) 
and consider the limit 



(0(x))^ = hm (0(x)) 



R 



lim 



\ a 


e^^i^O(0,0) e 


~Hl{R~-x) 


Bis) 


{Bk 


g—H^R 





(2.14) 



where Hl is the finite volume Hamiltonian and represent the boundary states in finite volume 

(see fig. 1(b)). It will be the subject of the next subsection to properly expand them in the basis of 
the eigenstates of Hi. 



2.1 Boundary states in finite volume 

The boundary states can be expanded in finite volume as 

where l^')^ are eigenstates of the finite volume Hamiltonian Hi. In particular the function G^{L) 
determines the large R behaviour of the partition function and can be written as 

Here fj is the contribution of a single boundary to the ground state energy Eq{R) and QjiL) is 
the standard non-perturbative ^-function"^ which was introduced in critical systems by Affleck and 
Ludwig [26]. In theories with only massive excitations in the bulk the ^-function decays exponen- 
tially. The exact integral series for the (7- function in massive theories was derived in [27]; for a first 
treatment of a non-trivial massless flow see [28]. 

The excited states |^')^ of a finite volume Hamiltonian can be described in a large volume as 
scattering states consisting of n particles with rapidities On given by the solution of the asymptotic 
Bethe-Yang equations 

Qk = niLsinhOk + ''^^S{6k — 0j) = 27rlk , k = l,...,n. (2.15) 

It is convenient to define a continuous two-particle phase shift function by^ 

S{e) = -e^^(^) , 6{-e) = -6{9) . 

*This should not be confused with the one-particle boundary coupling gj. 

^ We made use of the fact that the effective statistics is fermionic, i.e. ^(O) = —1. All known integrable models 
possess this property, the only counter-example being the free boson. In a theory with more than one particle species 
the most convenient choice for the definition of the phase shift is Sab{0) = Sab(0)e"^"'''®''. 
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This poses the following prescription for the momentum quantum numbers: 



/fc G Z for odd n , € Z + — for even n . 

The quantum numbers {/i, . . . , /„} completely determine the individual scattering states, therefore 
they may be used to label the sates as 

IUl---^}>L, 

where the states are normalised to unity: 

l({Ji . . . In}\{l'i . . . l'm})L = • • • Sj^j^ . 

The total energy and momentum are calculated additively as 

n n 

E = ^m cosh 9, + ©(e"^^) , p = ^ m sinh 9i + ©(e"^-^) . 

i=l 1=1 

The exponential corrections are governed by the mass scale /i which is uniform in the sense that it 
is determined by the analytic structure of the S-matrix and it does not depend on the particular 
multi-particle state in question. For a systematic treatment of exponential corrections to excitation 
energies see [29, 30] and also recent papers [31, 32]. In a massive field theory it is expected that the 
leading exponential corrections of more complicated quantities (like form factors [33] and correlation 
functions) are also of order e~^^. 

For later use we introduce the density of states (in rapidity-space) in the n-particle sector of the 
theory, which is given by the Jacobian 

Pn{9i,...,9n) = det Jkl , '^'^^ ^ ■ ^'^'^^^ 

Apart from normalisation issues the boundary states in finite and infinite volume should have 
the same structure. Therefore, it is natural to expand the finite volume boundary state up to the 
two-particle contribution as 

|s/> = GO(L)(|O)^ + |iVi(L)|{O})^ + 5]K,(0)iV2(0,L)|{-/,/})^ + ...) , (2.17) 

/>o 

where in the last term it is understood that the rapidities 9 are the solutions of the appropriate 
Bethe-Yang equation 

Qi{9) = mLsmh9 + S{29) = 2TrI . (2.18) 

Note also that the sum in (2.17) only runs over / > because the states |{/, —I})^ and |{— /, /})^ 
are identical. In (2.17) we introduced the normalisation factors Ni{L) and N2{9,L) which are not 
determined by first principles. In the following we calculate them including all finite size effects 
which scale as negative powers of L (only neglecting those which decay exponentially). To this 
order we consider vacuum expectation values in the i? — )• oo limit with fixed x, both in a large finite 
volume L and directly in the infinite system. 

In infinite volume the first three contributions read [13] 

(0(x)), = (O) + |iFOe-™^- + ^K*{9)F^{9, _0)e-2mxcoshe ^ _ _ _ _ (2.19) 
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The same quantity can also be evaluated in finite volume as 
(0(^))a = (0) + f A^iW({0}|0|0)^e-^o^ + j;/CWiV2(^,i^)({-/,/}|0|0)^e-^^- + ... , (2.20) 

where 

Eq = 771 + 0(6"^^) , Ei = 2mcoshe + C)(e-^^) , 

and the finite volume form factors are given by [15] 

The summation in (2.20) can be replaced by an integration leading to 

+ - [ ^Pl^^}^M^Kl{9)F^{e, _0)e-2™--he ^ 0(e-3™-) + ©(e"'^^) , (2.21) 

2 J 2-K ^p2{e,-e) 

where pi{9) is the constrained density of states given by 



dQ 

In a massive theory 



(e) = -^ = mL cosh 9 + 2ip{29) . 
d9 



{0{x))^-{0{x)t - Oie-^"^). 

Requiring that expressions (2.19) and (2.21) only differ by exponentially small terms for any scalar 
operator yields 



Ari(L) = ^^ + 0(6^^^^) , N2{9,L) = ^ '^f' +0(e-^^) . (2.22) 

Note that this argument essentially coincides with that used in [15] to determine the normalisation 
factors of the elementary finite volume form factors. 
The two-particle density satisfies 

P2{e,-e) = p^{9)-pi{e) , 

therefore 

We stress that the eqs. (2.22) include all finite volume corrections which behave as negative powers 
of L, in particular there can be no additional 0{1/L) corrections to N2- This will become important 
in the evaluation of vacuum expectation values at finite R. 

It is straightforward to extend the argument above to higher excited states. The finite volume 
boundary state is given up to the four-particle terms by 



\Bf) = G^j{L){\0)^ + ^-^V^m)^ + Y,Kjm2{e,L)\{I,-I})^ 

I 

+ E ^K,{9)N3{9, L)\{I, -1, 0})^ + ^Yl Kji9i)K,{92)N,iei,92, L)\{I, -I, J, -J}) A . (2.24) 



2 ' ^ ' ' '^'^ 2 
/ ij 
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The three-particle sector consists of states with rapidities {—^,0,^} where 9 is determined by the 
single quantisation condition 

03(6') = mLsinh6' + 5{e) + (5(26*) = 2t:I . (2.25) 
Therefore the normalisation of the three-particle states is given by 



where 

p^{e) = = mL cosh9 + ^{6) + 2(^(20) . 

do 

In the four-particle case the quantisation condition for the rapidities { — ^1, ^1, — ^2, ^2} is given 
by the system of equations 

Q4,i = mL sinh 01 + 5(^1 -^2) + 5(01 + ^2) + 5(2^1) = 2^/i , 
04,2 = mLsinh02 + 5(02 - 0i) + 5(0i + 62) + 5(202) = 27r/2 , 

yielding the normalisation 



iV4(01,02) = ^^^^^^^^^ 



^4(^1,02) 

where 

/O4(0i,02) = det J with Jik = , i,k = 1,2. 

2.2 Large volume expansion of the Casimir energy 

Before turning to the evaluation of the one-point function we consider the Casimir energy in the 
two-boundary setting. In [34] the exact non-perturbative value was found to be 

E^{R) = -\ I ^mcosh01og(l + i^:(0)if^(0)e-^(^)) , (2.26) 

where e(0) is the solution of the boundary TBA equation 

e{e) = 2mRcoshe- J ^99(0 - 0') log (l + K;(0')i^/3(0')e~'^^'^) , (2.27) 

where (p{0) = d5(0)/d0. Expression (2.26) is normalised to have the asymptotics Eq{R) — )• as 
i? — )• cx). The actual ground state energy also includes bulk and boundary contributions, which can 
be calculated from the boundary TBA [35]. 

Equations (2.26)-(2.27) were derived in [34] for the case = 9/3 = but it was argued that 
they yield the correct result even in the presence of zero-momentum particles. Although this claim 
was supported by several numerical checks [35], a rigorous proof is still missing. Moreover, it turned 
out that there are certain cases when (2.26)-(2.27) cannot be true. 

In [19, 20, 36] it was shown that if 7^ 7^ 5^ then the leading contribution in the integral in 
(2.26) is given by 

Eo{R) = -m^-MAe-^R + . . . . (2.28) 



10 



On the other hand, the authors argued that the leading behaviour of the Casimir energy is always 

Eo{R) = -m^e-""^ + ... , (2.29) 

irrespective of the sign of Qadp- This discrepancy shows that the validity of (2.26) is restricted to 
the case gaQp > 0. On the other hand, it was shown in [19, 20, 36] that it is possible to derive a 
proper analytic continuation to the region gaQ/s < with the correct leading behaviour (2.29). 

Here we derive a large R expansion of the Casimir energy by means of the finite volume regu- 
larisation of the partition function. Our results will be compared to (2.26). We follow closely the 
treatment of [19, 20], however we are able to extend their method to obtain the second order terms 
of the Casimir energy. In order to keep the exposition simple, we only consider the region Qadis > 
where the BTBA (2.26)-(2.27) is applicable. 

The partition function Z(L, R) can be evaluated in two different channels. Treating L as time 
variable and R as space one obtains in the large L limit 

Z = e--^o(«)^ (1 + Oie-""^)) , 

where Eq{R) is the ground-state energy in a finite volume R. It follows that 

Eo{R) = - hm . 

L— >oo Li 

Treating R as time variable (corresponding to thermal field theory with temperature T = 1/R) the 
partition function can be developed into a low-temperature expansion. In the regime mLe""^^ ^ 1 
it is sufficient to sum over the low-lying excitations: 

Z = l + N^Mle-^^R + (iV2(0))^i^*(e)K;3(^)e"^'"^™"^^ + Oie-^""^) . (2.30) 

Note that in the expression above we did not perform a complex conjugation for Qa, although the 
amplitude Ka{0) is conjugated. The motivation for this prescription is a delicate issue. In unitary 
theories the one-particle couplings are always real, therefore it is only the non-unitary case which 
has to be addressed. In non-unitary theories Qa is purely imaginary, therefore complex conjugation 
simply results in changing the sign of the appropriate term. On the other hand, the space of the 
states is not positive-definite, in fact one-particle states can be assigned a negative norm. This 
additional change of sign accounts for the correct prescription in (2.30). A different argument is 
provided by the corresponding conventions in CFT, where the prescription for the inner product 
does not include complex conjugation. Moreover, when comparing our form factor expansion with 
TCSA results in section 4.2 we will show that this choice of the sign is the correct one. 

If mL ^ 1 one can safely substitute the normalisations (2.22) into (2.30). The logarithm of the 
partition function is then written as 

= mM^e-^ + y ^^^K:(e)K,(0)e-2-^^-'^^ 

_ ^rn^L (M^e— ' + ©(e"^^) + ©(e'^-^) . (2.31) 

The first term is already in complete accordance with (2.29): in the following we also evaluate the 
second order terms. First of all note that the summation over / is divergent because of the double 
pole 
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In order to determine the leading singularity it suffices to use the lowest order approximation in 
(2.15) 

27r/ r _ 1 3 5 

2'2'2''" • 

The divergent part of the summation in (2.31) can then be expressed as 

.2 „2 / r N 2 „2 2 

-2mR 



32 



where we used the identity 



i\' f^V fbV 



Therefore the 0(Le~^™^) terms cancel exactly in (2.31), as needed to obtain a meaningful L — )• oo 
limit. Note that this result provides new evidence for the normalisation condition gj = gj/2, because 
the L — )• oo limit would not exist with a different choice of gj . 

In the following we also evaluate the finite left-over piece of the C)(e~2™^) terms in (2.31). In 
order to make calculations as transparent as possible we first introduce the following theorem: 

Theorem 1 Let 9j be the solutions of the quantisation condition 

Qi{e) = mL smhe + 5(26*) = 2ttI 

and let f{9) be a symmetric function which apart from a double pole at 9 = is analytic in a 
neighbourhood of the real axis: 

G 



Then the expression 



f{9)^-^ as 0^0. 



has a regular behaviour at large L with the L — )■ oo limit given by 



lim SiL) = If + Kf, 



where 



The proof can be extracted from the Appendix B of [37]. In addition, we also present a different 
derivation in Appendix A. In the present case Theorem 1 can be applied with the substitutions 



g2 2 

f{9) = COsh0K«(-0)i^^(0)e^2miJcoshe ^ ^ ^ ^^^-2mR 



Here we made use of the identity K^{9) = Ka{—9). It is easy to see that with this choice f{9) is 
indeed symmetric and it is analytic apart from the prescribed double pole at = 0. 
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The results of Theorem 1 for the Casimir energy yields 



2 2 

^ 9agi3 -mR -2mRVi^) 

Er)(R) = —m -e — m -e -^^'^"-J—^ 

4 4 4 



_ i y l^mcosh^ ^K„(-0)K/3(e)e-2"^^™^^^ - + 0(e-3-^) . (2.32) 

In order to compare this result to the BTBA equations it is necessary to derive a large R 
expansion of (2.26). First of all we regularise the integral along the lines of [19]: 

E(R) = ^ _ i / eosh S log ( ^^K„(-e)Kme-'<'A ,, 33, 

4 2 y 27r \ 1 I _E^EPp-e{0) j 

Using the same trick for the integral in (2.27) it is possible to derive the leading correction to the 
pseudo-energy: 

e{e) = 2mR cosh 6 - Ml^^^^y^raRcoshe ^ Q^^-'imR^ ^ (2.34) 

and in particular 

Ve--(O) = e-™^ ( 1 + Ml^(i^)e-^R + ...). (2.35) 



4 

The second term in (2.33) can be expanded as 

r,2 „2 



K„(-0)K«(0)e-^W - -^4|-:e-^(0) / 5^51 

4smh y 



(2.36) 



In the last step we used the fact that the extra term in the denominator only makes a difference for 
~ e"™^, therefore it only contributes to the higher order terms. In the last expression we may 
finally substitute the zeroth order approximation £{9) ~ 2mi2cosh0. Putting everything together 
we arrive at 



^2 „2 



_ i I ^mcosh^ |^i^,(-0)i^^(^)e-2-^^°^^^(^) - + 0(e-3-«) . (2.37) 

This is in complete accordance with (2.32) and can be regarded as a strong confirmation of the 
consistency of our calculations. Apart from checking the identity gj = gj/2 and also the validity of 
the BTBA equations, we also confirmed the normalisation factor N2{0,L). In particular, a different 
0{1/L) term in N2{9,L) would yield a different pre-factor for the term 93(0)6"^™^. 

Finally we note that our finite volume regularisation scheme is not sensitive to the sign of Qadp, 
in particular one arrives at the same result (2.37) also for gad/^ < 0. On the other hand, in this 
regime one has to use the modified BTBA equations given by eq. (2.7) of [20]. We leave it as an 
exercise to show that a careful large R expansion of this modified BTBA indeed reproduces (2.37). 
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2.3 Connection between boundary states and form factors 

As a remark to this section we would like to point out an interesting similarity between the elemen- 
tary finite volume form factors (0|O|^')^ and the amplitudes gf{L) defined by 

{B!f\^)^ = e-f^'^gfiL), j = a,/3, 

where fj is the contribution of a single boundary to the ground state energy on a strip. Both 
quantities measure the overlap of a normalised eigenstate \^)]^ with a non-normalisable extended 

state, namely the boundary state -Bj^^ and the state O\0)j^ created by acting with a local operator 
on the vacuum. To demonstrate the analogy we compare the main results in the case of the specific 
two-particle state l^')^, = |{/, — We have seen that 

gf{L) = K,{e)N2{e,L)g^{L) + ©(e"'^^) . 

The vacuum amplitude g^{L) only include exponentially small corrections (see below), therefore the 
excited state amplitude can be written as 

gJ{L) = K^{e)N2{e,L) + 0{e-^'^). 

On the other hand, the finite volume form factor was shown in [15] to be given by 

(0|O|^), = ^^^^^-^^ + 0(e-^^) . 

Vp2(^,-0) 

We observe that apart from the exponential corrections both objects are determined by the cor- 
responding infinite volume quantity and a normalisation factor which only depends on the finite 
volume density of states. In the following we make some comments about the possible structure of 
the exponential corrections. 

First of all it is instructive to recall the known results in the case of I'l')^, = lO)^^, which corre- 
sponds to the exact non-perturbative (7-function and to the finite volume vacuum expectation value 
{0)j^, respectively. In the case of the non-perturbative (/-function the final result was expressed in 
[27] as 

2\ogg%L) = \j ^{^m - m - 2^(20)) log (l + e'^^^)) 

where ^j{9) = —i-^log{Rj{d)) and e{9) is the solution of the periodic-boundary- conditions TBA 
equation 

e{6) = mL cosh 9- J ^^{9 - 9') log (l + e'^^^')) . 

In the case of the vacuum expectation value the relevant exact result is the LeClair-Mussardo series 
17) 

<°>r"" - E ^ / ^ ■ ■ ■ / It (ri ttW) 

n=0 \i=l / 



where e(^) is the same pseudo-energy function as above. There is a striking structural similarity 
between the two series, most importantly the weight functions appearing in the integrals are exactly 
the same. 
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Based on these similarities we conjecture that the structure of exponential corrections to the 
excited state quantities (O|0|^')j;^ and gf (L) will be similar to each other as well. It was remarked in 
[28] that the exact series for the amplitudes df {L) probably involves the solution of the correspond- 
ing excited state TBA. The appearance of the excited state TEA was also observed in preliminary 
studies of exponential corrections to form factors [38]. The exact result for gf {L) will probably 
involve an integral series similar to (2.38) normalised by a "dressed" form of N2{6,L). 



3 Evaluation of the vacuum expectation value 



In this section we develop a large R expansion of the vacuum expectation value (v.e.v.) {O)'^^ . 
We use the finite volume regularisation scheme described in the previous section and express the 
expectation value as 





^-Hlxq ^-Hl{R-x) 


^1) 


{Bk 







(3.1) 



We will calculate explicitly the v.e.v. with the boundary states truncated to contributions with up 
to four particles. In order to simplify notations we define 



where 



Br, 








B^ 


)« = 






B]) 


)(^) = 


bT. 













Lb 



|0)i, 



n=0 



j = a, /3 



|iVi(L)|{0})^ , 
j;K,(^)iV2(0,L)|{/,-/})^, 

j;|A-,(^)iV3(0,L)|{/,-/,O}), 

\ K,{er)K,{e2)Ni{er, e^^mi, -i, j, - j}), 



(3.2) 



IJ 



We also define 



Cnm = (")(i?^|e-^^-0(x)e- 



~HL(R~x)\^L\i^) 



13, 



In the following we perform a double expansion of (3.1) with expansion parameters e~"^^ and 
^-m(R-x) jj-^ order to keep track of the various terms let us introduce two auxiliary variables u and 
V. They are use to count orders of e~"*^ and g-'^C^-^) and at the end of the calculations both are 
set to 1. Then the v.e.v. takes the form 



We define 



(3.3) 
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where the first terms are given by 

Zo = 1 , Zi = M^mLe-'"^ , (3.4) 

and 

Z2 = Y^ e-^'"^°''^^^K*^{e)Ki^{e)N2{e, if . (3.5) 

I 

The inverse of the partition function is expanded as 



where the first few terms read 



■^0 — 1 ) Zi — —Z\ , Z2 — Z\ — Z2 



Putting this together we obtain 



i>nm — ^ Cn-l,m-lZl , (3.6) 



, . . . , 



I 

where the first few non-trivial terms are given by 

Dim = Cim — ZiCo^m-l , m = 1,2 

D2m = C2m — ZiCi^m-l + {Zf — Z2)Co,m.-2 , m = 2, 3, . . . . 

We expect that the quantities -Dnm have a regular behaviour as L — t- oo and for the actual limit 
we define 

Dnm = lim Dnm ■ (3.7) 



The vacuum expectation value is then expressed as 



mm 

n,m, 



Note that the individual terms contributing to Dnm may contain divergent pieces which scale with 
positive powers of L. The most singular terms in the n particle sector carry a factor of (mLe"*"^)", 
therefore the expansion is valid in the regime 



The L — )■ oo limit of the complete series (3.3) is to be understood as an analytic continuation. 

The evaluation of the individual Cnm is built on the knowledge of the finite volume form factors. 
In the general case they are given by [15] 

({/„..., In}\0\{Ju Jm}), = + , (3.8) 

where it is understood that the rapidities {^i, . . . , and {9'i, . . . , O'^^} are solutions to the corre- 
sponding Bethe-Yang equations. Formula (3.8) is valid whenever there are no coinciding rapidities 
in the two scattering states. In [15] it was shown that the only two situations when coinciding ra- 
pidities occur are the case of diagonal form factors and matrix elements between parity-symmetric 
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states containing zero-momentum particles. The general rule for evaluating disconnected contribu- 
tions can be found in [11]; in the present work we will cite the necessary results whenever they are 
needed. 

Before turning to the actual calculations we have to address a very important question concerning 
the phase of the form factors. Note that equation (3.8) is to be understood up to a phase factor; in 
particular the order of the rapidities is not determined by first principles. In fact, an additional phase 
factor corresponds to a redefinition of the basis vectors, and this phase drops out from the calculation 
of some physical quantities, for example correlation functions in infinite volume. However, phase 
factors are utterly relevant in the present case because they do affect the final result for the v.e.v. 
A guideline can be established from the known results in infinite volume, where the correct phase 
is fixed by consistency arguments. First of all, the two-particle contributions to the boundary state 
are written as 

K,{e)A{-e)A{e) . 

This expression is symmetric in 9 as required by consistency. When it comes to the evaluation of 
the v.e.v. in the presence of a single boundary, the above definition yields the expressions 

...Kp{e)F^{-e,e,...) and ...Ki{e)F^{e + m,-e + i^,...) , 

where the dots stand for possible additional rapidities and amplitudes. Observe that both expres- 
sions above are symmetric in 6, therefore we apply the following rule: whenever there appear the 
amplitudes Kp{9i) and K^{9i) with some 0i, the explicit form and order of the rapidities substituted 
into the relevant form factor is given by {—Oi, 9i) and (0,, -l-ivr, —9i + respectively. The exchange 
of two pairs of rapidities does not make a difference, therefore the phase of the form factors is 
completely fixed by the above rule. 

The complex conjugation of the amplitude can be avoided by means of the identity 

Kl{9i) = K^{-9i) . 

An additional change of variables 0j — ?• —9i then amounts to the convention 

...K^{9)F^{-9 + m,+9 + m,...) . 

Note also that the presence of zero-momentum particles does not produce any additional ambiguities, 
because it is only the expressions 

K,i9)F^i-9, 9,0,...) = K,{9)F^iO, -9,9,...) 

which are symmetric in 9, i.e. the zero-momentum particle cannot be placed between the particles 
A{9) and A{-9). 

In the following we evaluate all Dnm with n + m < 4. We put forward that all our results are 
collected in Appendix C together with a pictorial representation of the individual terms. 

The contributions Dno and Dom with n,m = 1 ... 4 only depend on one of the boundaries and 
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they are the same as obtained in [13]: 

n _ rpO„-mx 



2 

o 7 ZTT zvr 

i)04 = ^/ ^^i^/.(0l)A;,(02)i^f (-01, 01,-02, 02)e-^™(-^^^^+^°^^^^) . 

Note that these integrals remain well-defined even if there are poles in the amplitudes Kj{9), because 
the form factors possess the appropriate number of zeros at 0j = as a consequence of the exchange 
axiom (2.9). 

The first contribution to contain a divergent piece is Dn which is given by 



where 



Dii = hm Cii-ZiCoo , (3.9) 

Cii = ^mL ({0}|O|{0})^ e— « . (3.10) 
The diagonal one-particle form factor for a generic / is given by [11] 

({7}ioi{j}>,= ^°'",;;''" +(o). (3.11) 

The two-particle form factor appearing in the expression above is free of divergences and by Lorentz- 
invariance it does not depend on 0. In the following it will be denoted by -F^^g. Specifying (3.11) to 
1 = one finds 

Dn = hm (Cn - Z^Coo) = M^F^^ e^^^ . (3.12) 

3.1 Evaluation of D21 and D12 

We first consider 

D21 = lim (C21 - ZiCio) (3.13) 

with 

C21 = f ^aZ^±!^li±^e-^--(^-)K.(0) , (3.14) 

where it is understood that the rapidities are solutions of the corresponding Bethel ang equations 

mLsinh0 + 5(0) + 5(20) = 2ttI , 
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and the constrained density is 

pi(6») = mLcosh6' + 299(26*) . 

Note also that we have aheady made use of the relations (2.22). Naively one would take the L — )• 00 
limit simply by replacing the summation with and integration leading to 

2 2 J 2tt 

However, the expression above is ill-defined at 9 = due to the poles of F3 and Ka- Not surprisingly, 
this divergence gets cancelled by Z\C\q as we will show shortly. 

It follows from the form factor axioms (2.9) and (2.11) that near Q\ 2 = 

F^(Q^ + ivr, 02 + ivr, 0) ^ 2i - 1^ (3.16) 



and therefore 

F^{-e + i^,e + m,Q)^^i F^\ + ... . (3.17) 

u 

The singular behaviour of the sum (3.14) is determined by the states near = 0. We use the 
approximation of Section 2.2 to determine the leading divergence as 

. P? {9a? 913 -ra(x+R) V /^l V = 4^ {9cd'^9p ( -rn{x+R) ^ /o i 

mL 4 V 2vr y j \I J mL 4 \2ttJ 2' ^' ^ 

This coincides with ZiCiq therefore (3.13) has a finite L — )• oo limit indeed. 

Similarly to the calculations of Section 2.2 we use Theorem 1 to obtain the finite-left over piece 
in (3.13). We apply the substitutions 

f{0) = F3^(-0 + i7r,0 + i7r,O)e-2'"™«hto-m{i?-x)^^(^) ^ G = 2{g^fF^e-"'^''+^^ . 

It is easy to see that f{9) is indeed a symmetric function, therefore its only singularity near the real 
axis is a double pole at 9 = 0. The net result is then expressed as 

il f ^ f fO(_q,- Vi\K (Q\p-2mcoshex-m(R-x) _ '^{9a? COsh 9 (j^^^ 



D2i = ^ —I fO{-9 + m,9 + ivr, o)ir„(0)e-2™--^''"-"^(^-") - "^^"^ M^^^^^-^ ^-^fl+x) 
4 J 27r \ sinh 9 

+ e-(-+^)ff,(5.)2Ff^. (3.19) 
A similar calculation with the roles of the two boundaries exchanged yields 

9^ f ^ f pOf- a n\T^ ^-2771 cosh e{R-x)-mx _ ^ (g/3 ) ^ COsh 6* (^fl- 



D,2 = ^ —i F3«(i7r,-0,0)i^«(0)e-2™-^he(ii-x)-mx _ ^UWM^:^g-m(2fl-x) 
4 J 27r V sinh 9 



+ e-^'^-^^gMpfF?^. (3.20) 



3.2 Evaluation of D13 and D31 

We first consider D^i = limi_j.oo(C'3i — ■Z^iC'20) where 



C31 - = Y,^3{0)V^^K{9) ({/,-/, 0}|O|{0})^ ^-mR-2moosh9 x_ 



4 

9a9p 



mLe— ^5]iV2(0)({J,-J}|O|O)^e-^^-K*(0) 
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The three-particle normahsation is given by 

Vp3(^,-^,o) 



The four-particle matrix element above includes a disconnected term due to the zero-momentum 
particles which is given by [11] 

(U -/ OllOliO)) = + -0 + 0) + mLFOje, -6) 

^/p-i{e,-e,Q)mL 

while the connected part of C31 can be transformed into the integral 

9a9l^-mR I ^j^^^g^jpO^_0 ^i^^g ^i^j^^Q^^-2mcoshe X _ (3 21) 

8 J 211 

It follows from the form factor axioms that F^[—6 + iTT,9 + iTT,iTT,0) is well-defined [11] and has a 
regular behaviour as a function of 9. Therefore we can apply the exchange axiom (2.9) to show that 

lim F9(-e + iTr,e + ivr, ivr, 0) = , (3.22) 

(9-5.0 

therefore the integral (3.21) is regular even at = 0. 
The divergent pieces of D^i are given by 

Note the summations over / and J run over the solutions of the different quantisation conditions 
(2.18) and (2.25), respectively. However, the density factors are the corresponding ones in both 
sum, therefore it is allowed to replace the summation with an integration in both terms. Doing this 
we make errors of order e~^^ because both integrands are regular, therefore the above expression 
vanishes in the L — )• oo limit and the net result is 

= M^e-™« j ^KUe)F^i-e + zn,e + ivr, ivr, 0)e"2™h. . ^ (3 23) 

and similarly, 

D^s = Mle-^R f ^i^.(0)Ff (-0 + i7r,0 + m,i7r,O)e-2™-«he (R-x) ^ 
8 J 271 

3.3 Evaluation of D22 

Here we consider 

D22 = lim (C22 - Z2 (O) - ZiCii + {Zif (O) ) . (3.24) 
Evaluating C22 one has to treat the diagonal and off-diagonal matrix elements separately: 

C22 = J2N2{ei)N2{e2){{i,-i}\o\{j,-j})Le-'''^-''^<^''-^'>K^{-ei)K^{e2)+ 

+ Y.N2{ef{{i,-i}\o\{i,-i})^K^{-e)Kp{e)e-'''^ . (3.25) 
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Here it is understood that the rapidities 6i and 62 solve the quantisation condition (2.18) with 
momentum quantum numbers / and J, respectively (including the diagonal case with 9i = 62 and 
I = J). In the off-diagonal case one has 

N2{e,)N2{e2){{i, -i}\o\{j, -J}), = + + _ 

It would be desirable to convert the sum over these matrix elements into an integral. However, one 
has to be careful because of the ambiguity of the form factors near 61 = 62 and the possible poles 
of the integrand at = ^2 = 0. As a first step we write 



E -^f (^1 + iT^, —Ol + ^TT, —02, O2) , a \ (a N„-2m cosh 6»ix-2m cosh 6l2(-R-x) 



]_ jdBl j (^^)^0(^^+ .^^.^^^ .^^_^^^^^)g-2rncoshe,x-2™cosh^?.(iJ-x) 

_ ^ Ff(0 + .vr,-^ + ^vr,-0,0) ^^^_^^^^^^^^_,^^ ^ ^^^^^^ 

where F^{0 + in, —6 + ivr, —9, 6) is the continuation of the function F^{6i + ivr, —61 + m, —62,02) 
io 9i = 62 = . Kt first sight there seems to be a double kinematic pole, however this is not quite 
true. In [11] it was shown that diagonal form factors may be assigned a finite value, although the 
diagonal limit will depend on the particular evaluation method used. In the following we restate the 
results of [11] which are needed to evaluate (3.27). 
Consider the quantity 

FZi^i + ei + ivr, . . . , 0„ + e„ + ivr, 6'„,, . . . , 0i) , (3.28) 

where the singularities have been shifted off by the infinitesimal quantities e^. It was proven in [11] 
that there exists a finite limit when all go to zero simultaneously with their ratios fixed. Moreover, 
there are two special evaluation schemes which respect the physical requirement that diagonal form 
factors should not depend on the order of the rapidities. First of all, one can consider the symmetric 
limit with £1 = £2 = ■ ■ ■ = e: 

Fg^ ,(01, . . . , 0„) = lim FgXOi + e + ivr, . . . , 0„ + e + ivr, 0„, . . . , 0i) . (3.29) 

On the other hand, one can also consider the connected part of the diagonal form factor which is 
defined to be the contribution to (3.28) which does not contain any singular factors of the form 
£i/£j and products thereof: 

-^2n,c(^i> ■■■,0n) = (finite part of) ^2^(01 + ei + ivr, . . . , 6l„ + e„ + ivr, 6l„, . . . , ^i) . 

The general structure of the singularities in (3.28) together with the relation between the sym- 
metric and connected evaluation schemes was worked out in [11]. Here we only need the formula 
relevant to the two-particle states, which states that given two infinitesimal numbers £1 and £2 one 
has 

F2{ei + ei + ivr, -62 + £2 + ivr, -02, ^i) = ifc(ei, ^2) + {- + -\ ^{Oi - 02)if . • (3.30) 
In particular, the symmetric evaluation with £1 = £2 is given by 

F« (01, 02) = f2A0i, O2) + 2(^(01 - 02)^2^, . (3.31) 
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In the present case we need to evaluate F^{6i + iir^—Oi + ivr,— ^25^2) near 61 = 62 = which 
corresponds to an "antisymmetric evaluation" with £2 = — ei: 

F^{e + ivr, -e + ivr, -9, 6) = lim Ff (0 + e + ivr, -0 - e + ivr, -0, 0) 

(3.32) 

= FO(e,-0)-299(2e)F2« . 

Note that the integral in (3.27) is well-defined even if there are poles in the factors Ka[—0)Kp{0) 
because ^4^(6*1 + ztt, — 6*1 + itt, —62,02) possesses a double zero at ^1 = ^2 = (for a proof see 
Appendix B). 

In order to evaluate C22 we also need the diagonal finite volume form factors. A generic diagonal 
four-particle matrix element reads [11] 

F4M,02) + (pi(9i) + Pi{B2))F?, 
{{h,l2]\0\{hj2])L= ' Z L ^ +{0) . (3.33) 

P2(,fc'l, t72j 

Specifying the above formula to the present case one obtains 

({-/,/}|0|{-/,/})^ = ' ' + (O) . (3.34) 

Substituting the results of (3.27), (3.32) and (3.34) into (3.25) we obtain 
(c22-Z2{0)-ZiCii + {Zif{0) 



^ ^K^i-(^i)K^i92)F^{ei + i7T,-ei+ in, -02,02)6 



2m cosh 0ix—2'm cosh 82{R—x) 



+ F2Z [Y. -r^fo.{-e)Kp{e)e-^^^^^'^<^^ - i^^^mLe-2-«j + 0{e-^^) , (3.35) 

where we applied the normalisation given by (2.22) and we also made use of the relation (3.31) and 
the identity 

pi{0) = p{0) + 2ip{20) . (3.36) 

The summation in the third line of (3.35) still contains a divergent 0{L) piece. However, it can 
be proven along the lines of the previous subsections, that this divergence gets exactly cancelled. 
One can use Theorem 1 with the substitutions 

f{0) = 2A',(-^)K;3(^)e-2-™^h9R ^ Q ^ i9a^^.2mR 

to obtain the finite left-over piece. The net result reads 

I?22 = lim (C22 - Z2 {O) - ZiCn + (Zi)' (O) ) = /22 + J22 , 

where 

/22 = 7 / -^-^i^«(0i)i^/3(02)i^f (-^1 +ivr,0i +i7r,-02,^2)e-2™--^^^"-2"^™^'^^^(^-") (3.37) 
4 J zTT zvr 

and 

^22 = ^2^. / '!^{K^{-9)K,{9)e-^--^^'^ -^^^fl^ 

' J ZTT \ 4 smh U / ' o 

(3.38) 

With this we have finished the calculation of the terms in (3.3) with n + m < A. The reader can 
find the various terms collected in Appendix C together with a pictorial representation of them. 
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3.4 Connections to other problems 



In this subsection we compare our results to those obtained in previous works. Also, we establish a 
connection to the problem of finite temperature correlation functions. 



3.4.1 Comparison with a proposal of LeClair et al. 

The first result in the literature to deal with expectation values on a strip appeared in [4]. In 
Appendix C the authors conjecture the general result to be 

oo ^ 

iOi^))R = EE VAL^- ' (3.39) 



n=0 £i=± 



where 



I^Le. = n {/^ S^^»(^^')} ^2n(^l' -^1' • • • ' -On)e„...,e. • (3.40) 

The weight functions above are defined by 

„— 2mxcoshS TV' / n\ —2m(R—x) coshO ry- ( q\ 

and the form factors are given by 

-^2n(^l> -^1) • • • -0n)sl,...,er^ = F^ni^l + ^^1' "^1 + ^^1) • • • , 6'n + i^n, -On + i^n) 

with ai = 7r(l + ei)/2. It is assumed that Qa = Qp = 0. 

It is easy to see that some of our results are reproduced by the above series. However, there is 
also a serious discrepancy: (3.39) does not include the diagonal contribution J22 given by (3.38). 
This can be proven by putting x = R/2 and considering the large R limit of (3.39). In this case the 
nth term behaves as 0{e~^^^) while J22 is of order e"^*"^. Therefore J22 can only appear in the 
terms n < 2. On the other hand these contributions yield 

J« = D20 + ©(e"^'"^) , = D02 + 0(6-3"^^) , 

= + 0(6-4"^^) , /i'l = Do4 + 0(e-4"^«) , 
= j(2| ^ j^^ + 0(e-4-«) . 

All the above terms are well-defined and finite, thus it follows that J22 is not included in (3.39). 

We will demonstrate in section 4.2 with the help of non-perturbative numerical results that the 
term J22 must be included in the expansion, thus the series (3.39) can not be complete in the general 
case. 



3.4.2 Relation to quench problems 

In recent papers [2, 3] the integrable QFT approach was applied to quench problems. In particular, 
they considered the time evolution of expectation values of local operators after a sudden global 
quench which changes the Hamiltonian from Hq to H, where H (possibly also Hq) is considered to 
be integrable. The main assumption is that the initial state of the system (which is the ground state 
of Hq) can be expanded in the multi-particle basis of the integrable Hamiltonian H as a boundary 
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state. In the simplest case with only one particle type in the spectrum the corresponding expression 
is 

|i?) = exp (I ^K{0)A{-9)Ai9)^ |0) , 

where K(9) is an arbitrary function satisfying K{0) = S{26)K{—9). The time evolution of an 
expectation value is then given by 

(0(0, t)) = {B\e'"'0{Q,Q)e^'^'\B) . 

It is straightforward to establish a connection between the real-time evolution of the expectation 
value (0(0, t)) and the static (Euclidean) quantity {0{x))^ considered in this work. In fact, the 
quench problem above corresponds to the choice 

K^{9) = Kp{9) = K{9) , 

and the expectation value (0(0, t)) is obtained by the analytic continuation 

R = 2to , X = To — it , 

where tq is a constant which is analogous to the extrapolation length introduced in the field theory 
approach to boundary critical phenomena [1, 3]^. 

The results of [2] can be reproduced by our expansion with the analytic continuation explained 
above. The main difference between [2] and the present work is that [2] concerns the sine-Gordon 
model which is a non-diagonal scattering theory. However, the authors calculate contributions only 
from low-lying breather states which possess diagonal scattering. Therefore our methods directly 
apply to the situation in [2]. Moreover, our expansion includes higher order terms which were not 
considered in [2]. 

In the following we also compare our results to those of [3], which considers the i — )• oo limit of 
the expectation value (0(0, t)). For this quantity the authors obtain the series expansion 

(o) - 1 ^ n {/ f twf} ^&=<«" • • • • • 

where the weight function is given by 

0(61) = e-2'"^«™^''^i^(0) . 

The key idea of [3] is that in the t — t- oo limit only those contributions survive which do not depend 
on t. In our Euclidean setting this means, that the relevant terms do not depend on x. If we were 
able to derive the complete spectral expansion for {0{x))'^ , then in principle it would be possible 
to reproduce (3.42) by collecting those terms which are completely independent of x. However, we 
can already make a definite statement about the n = 1 term in (3.42). It is easy to see that up to 
higher order contributions this term is equal to our J22 given by (3.38), which is indeed independent 
of X. This way we have established the first rigorous support for the series (3.42)'^'. However, at 
present we can not make any statements about the higher order terms. In particular only a detailed 
study of the higher order terms in {0{x))'^ can decide whether some kind of a "dressing" of G{9) 
is needed to obtain the correct result, as it was suggested in [3]. 

* Keeping a finite (non-zero) tq is necessary to normalise the boundary states properly, otherwise the spectral series 
will not be convergent. 

^From a different point of view, the presence of the term J22 in (3.42) supports our claim that the series (3.39) 
(which does not include J22) can not be complete. 
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3.4.3 Relation to finite temperature correlation functions 

There is a way to connect the vacuum expectation values considered in this work to the evaluation of 
thermal correlation functions. The relation is not physical in the sense that it can be established only 
on a formal level. However, it is certainly worthwhile to explore this correspondence. For simplicity 
we suppose Qa = 913 = 0. Let us define a (non-local) operator B by specifying its form factors: 

{9i,92, . . . , 92n-l,92n\B\d'i,62, . . . , ^2m-li ^2m) = 

= (6{ei + 02)i^a(-ei)) . . . (<5(02n-l + e2n)i^a (-^2n-l )) 

X [6{e[ + e',)Kp{e[)) . . . [6{e',^_, + 9',jK^{e',^^,)) + (•.•). 

where the dots in parentheses stand for disconnected terms which appear if two pairs (02i-i5^2i) 
and (^2j-i'^2j) coincide for some i = 1, . . . ,n and j = 1, . . . ,m. The form factors with an odd 
number of particles on either sides are identically zero. Note that this definition of B incorporates 
information about both boundaries; actually the operator itself can be visualized by gluing the two 
boundaries together, thus forming a cylinder. 

With this definition of B the vacuum expectation value {0{x))'^ can be expressed on a formal 
level as a two-point function at finite temperature T = 1 /R: 

{0{x))l^ = (fi(0)O(x))r°'^^ . 

This correspondence implies, that higher order terms in the v.e.v. {Dij with i+ j > 4) will present 
technical difficulties which are very similar to those encountered in the evaluation of the thermal 
two-point function [7, 9, 17]. However, the exact relation between the two objects can only be 
established if the disconnected terms in the form factors of B are properly defined^, which is out of 
the scope of the present work. 

4 Comparison with TCSA 

In this section we compare for various cases our spectral expansion (3.3) computed in the previous 
section with the numerical results of the truncated conformal space approach (TCSA). We choose 
the scaling Lee^ang model, because despite being one of the simplest integrable models it already 
features many properties that we have discussed so far: it contains a single massive particle, it 
has tunable integrable boundary conditions with different gj values and all the form factors of its 
elementary field are known. Moreover, the TCSA is quite convergent for this model, which together 
with the properties listed above makes it an ideal testing ground for our method. In the first 
subsection we define the model and collect the relevant formulae which we will use to make the 
comparison in the second subsection. 

4.1 The boundary scaling Lee^ang model 

In this section we summarise the relevant properties of the scaling Lee^ang model described as a 
perturbed conformal field theory. We follow closely the treatment and notations of [13]. 

naive substitution of the form factors of B without any disconnected terms results in a series similar to (3.39). 
This way one misses the diagonal contributions like J22. 
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4.1.1 The critical Lee Yang model as a minimal conformal field theory 

The Lee^ang model [39] is the simplest example of a non- unitary conformal field theory. It corre- 
sponds to the M2,5 minimal model with central charge —22/5 and effective central charge 2/5. It 
has only two highest weight representations of the Virasoro algebra of weights and —1/5, so there 
are only two bulk primary fields, the identity IL of weight and (/? of weight —2/5. 

The conformal bounary conditions of the model were classified by Cardy in [40] and the field 
content on each of these is given by solving the consistency conditions written in [41, 42]. The 
model has only two conformally invariant boundary conditions, denoted by 1 and There is only 
one non-trivial boundary operator (/> = (p^^y^ that can live on the <I> boundary and there are two 

boundary changing operators, Tp = (A!^i*5 and tp^ = (p^^y^ , which interpolate between the boundary 
conditions 11 and These three boundary fields all have weight = = 1/5. 

For the sake of completeness we list all the operator product expansions (OPEs) and structure 
constants of the theory. The bulk OPE is 



the boundary OPEs are 
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and the two bulk-boundary OPEs read 

^{z)W = ''X\2{z-w) 
^iz)\^ = (*)i?^|2(z-ti;) 

The structure constants can be chosen as [43] 
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(4.1) 



We will need the conformal expectation value of the bulk field ip across the strip of width R 
with various combinations of boundary conditions: (1,1), (<!*, 1) and (<I>,<^). Our coordinates are 
< X < R across the strip and y running along the strip, and we normalise all our correlation 
functions so that the expectation value of the identity operator is always one. 

The correlation functions on the strip with boundary conditions a and /3 can be determined by 
mapping the strip to the unit disc or the upper half plane and inserting the appropriate boundary 



fields 



Under the exponential map z = exp(i7r(x — iy)/R the strip can be mapped to the 
upper half plane where the boundaries are mapped onto the negative and positive real axis. The 



-1/5- 
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boundaries meet at the origin (and at infinity) where boundary operators are inserted according to 
the boundary conditions. Through radial quantisation the Hilbert space of the theory is built on 
the representations corresponding to the boundary operators at the origin^. It follows that for the 
(1,1) boundaries the Hilbert space consists of the 1 module, in the (1,<I>) case of the ^ module 
and for (<I>, $) both modules are present with the ground state corresponding to (p. 

If at least one of the boundaries is the $ then the boundary field 0(y) may exist on it. For these 
cases the boundary one-point functions are 

(Hy)),^.) = (^) , (</'(y)>(. = f?! c^/ ■ (4-2) 



The one-point functions of the bulk field ip{x) on the strip correspond in general to chiral 
three-point functions on the upper half plane, and by the doubling trick they can be thought of as 
four-point functions on the full complex plane. Since the representation (j) has a null- vector at the 
second level, the correlation functions satisfy certain differential equations, the solutions of which 
can be expressed in terms of the four strip chiral block functions 

=(^)''Vi(|l,|j;fli;-tan^^), (4.3) 
h{d) = (2sin^?)2/5 , 
U{^) = (2sin^?)^/^ , 

where •& = xtt/R. The particular solutions are fixed by the boundary conditions, that is the corre- 
sponding bulk-boundary OPEs: 

/ p\ 2/5 



l<^'K/i(^), x> 



The latter is the unique combination of the chiral blocks with the proper asymptotic behaviour near 
the boundaries, for 1? — t- and — > vr. 

Since for the ($, $) boundary conditions both representations are present in the Hilbert space, 
we will need also the non-diagonal matrix elements 



{^Hx,y)\l) = '"KM^)= '^K{^^^^^) ' (4-5) 



Due to translational invariance in the "time" direction none of the expectation values depends on y. 



®Thus the states on the strip, unhke the bulk case, are characterised by their behaviour under one copy of the 
Virasoro algebra. 



27 



4.1.2 The scaling Lee^ang model 

The scahng Lee^ang (SLY) model can be defined as a bulk perturbation of the critical Lee^ang 
model by the term 

X dx dy(p{x,y) . (4.7) 



The resulting theory is an integrable massive scattering theory, containing a single particle with 
two-particle S-matrix [44] 

s{e) = -(i)(2) , (x) = :)i^:J, . (4.8) 

The mass M of the particle is related to the strength A of the perturbation by [45, 46] 

M = kX'/'' , K = 219/12^ = 2.642944 .... (4.9) 

^ 55/i6r(2/3)r(5/6) ^ ' 

The form factors of the bulk field ip were first computed in [47, 48], but following [13] we use 
the conventions of [49], with the difference that for us ^{x) is a real field. The function Fn can be 
parametrised as 

Fn{9l,...,en) = HnQnixi,...,Xn) fj H^^LllA ^ (4.10) 

. ^ . -^i \ -^7 

l<j 

where Xi = e^K The various terms m (4.10) are 

cosn n — 1 

^^^^ = cosh g + 1/2 ''^'^ " "'^"^'^ + ' ^'^■^^^ 



where v is given by 



10 



N 



= n 



n=l 



(^ + n + l/2)(^+n-l/6)(^+n-l/3)' 



(^ + n-l/2)(^+n + l/6)(^+n + l/3)_ 

^ / I- sinh(t/2)sinh(t/3)sinh(t/6) ^ ^ _ N (4.12) 

V JO tsinh^(t) / 

v{^) = 1.111544045... , 

and the normalisation factor is [50] 

Hn = ^[ 7= . (4.13) 

Here is the expectation value (ip) in the bulk: 

9 , 1 . 36 2 

-3l0r(i^M-5 ^ .9/^ . . 

Ip = ji-^^^ = (-1.239394325... Uf-^/^ . 4.14) 

(27r)-54r(i)r(|) 



■"This form with the finite product in front improves the convergence of the integral [49]. 
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The symmetric polynomials Qnixi, ■ ■ ■ have degree n(n — l)/2 and partial degree n — 1. They 
can be neatly expressed as a determinant of a matrix in symmetric polynomials [47, 48, 49]. Here 
we will need only the explicit form of the first few: 

Qo = 1 , Qi = 1 , Q2 = , Qa = 4'^^'^ , = <^?<^2^^i^ , (4-15) 

(n) 

where the elementary symmetric polynomials in n variables af- are defined by 

n n 

n(i+px.)=5^At" 



k=Q 



Let us turn now to the boundary conditions for the SLY model which were analysed in de- 
tail in [35]. The integrable boundary conditions are the 1 conformal boundary condition and the 
perturbation of the conformal <l> boundary by 



h J (j){x) dx (4.16) 
which we will denote by ^(h). The exact reflection factors for these two boundary conditions are 

R^^h){0) = Rb{e) , R^{e) = Ro{e), (4.17) 

where 

Rm = (i) (1) (I)"' {S{e + i7r!^)S{9 - • (4.18) 

The relation between b and h is given by [35] 

h{b) = - \hc\ sin(7r(6 + l/2)/5) . (4.19) 

For numerical calculations it is more suitable to use instead the dimensionless quantity 

which was determined in [51]: 

sin 2^ /rf2)\^/^ 

/ierit = -^'^'^'^'^'^'^j^^fJT^ ( f^j = -0.68528998399118... . (4.20) 

Finally, let us consider the boundary-particle couplings for the various boundary conditions. 
We recall that these quantities are defined via the residue of the reflection factor Ra{0) at the 
rapidity 6 = m 12 

^ ' 2 e-m/2 ^ ' 

In the SLY model we have 

51 = -i 2 \J 2/T- 3 (4.22) 

for the 1 boundary, and 

tan((5 + 2)7r/12) 
^^^^^= tan((6-2)vr/12) 

for the ^{h{b)) boundary. Notice that although the reflection factors for the boundaries l and 
$(/i(0)) are identical, the corresponding boundary-particle couplings differ by a sign. This means 
that the boundary states for these two boundary conditions differ only in the sign of the contributions 
from states of odd particle number. 
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4.2 The comparison 

4.2.1 Expectation values in the Truncated Conformal Space Approach 

The TCSA method aUows for calculating the spectrum and other physical quantities of a perturbed 
conformal field theory. The approach was developed in [52] for the bulk scaling Lee^ang model 
itself and later it was generalised to systems with boundary [35, 51]. The idea of the method is 
very simple: using the conformal field theory techniques one can calculate in the conformal basis 
all the matrix elements of the perturbing operator(s) and eventually of the perturbed Hamiltonian. 
If the conformal Hilbert space is truncated at some energy (or at some conformal level), the space 
of states becomes finite dimensional and the calculation of the spectrum and of the eigenvectors 
of the perturbed Hamiltonian boils down to the diagonalisation of a finite numerical matrix. This 
approach can be regarded as a close relative of the standard variational method, and although the 
errors caused by the truncation can not be easily controlled, it is generally believed, and checked of 
course, that they decrease with increasing the cut. We refer the reader interested in the details of 
the behaviour of TCSA to [53, 54]. 

The Hamiltonian of the boundary scaling LeeAfang model reads 

H{R; X,hi,hr) = Ho + X [ ip{x, y) dx + hi(j){x=0) + hrCpix^R) , (4.24) 

Jo 

where Ho is the conformal Hamiltonian and we have allowed for boundary perturbations at the edges 
of the strip (they can be present only for ^ boundaries). The matrix elements of the Hamiltonian, 
calculated by mapping the strip on the upper half plane are given by 



TT / C 



+ O''^' {G~'B),, + XI {^y^' {G-'Bi\j + xr {^y^' {G-^Br).,) , (4.25) 



where 



{Bi),, = mim , m^j = mi-m , (4.26) 

(i?),, = (i| r d^^ie'^^U , (4.27) 

JO 

and the conformal metric G defined as Gij = is needed because the basis vectors are not 

orthonormal. Here every operator is understood to be on the upper half plane, •& = ttx/R and we 
expressed everything in dimensionless form: 

r = MR , k! = K^2/5 = 0.097048456298 . . . , (4.28) 
XI = sin(^(6; + l/2)/5) /icrit , Xr = sin(7r(6^ + l/2)/5) /icrit , (4.29) 

implying that the energy eigenvalues are also measured in the units of M. 

Once the spectrum is known the expectation value of 99 can be estimated by 

where ^ = Mx is the dimensionless position of the operator on the strip and | $7 ) is the ground 
state eigenvector. If we know the matrix elements of ip{'&) on the upper half plane, the calculation 
of the v.e.v. amounts to matrix-vector multiplication. 

An efficient way to calculate the matrix elements is given in the Appendix of [55]. 
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(a) ($(-2), $(-2)) at r = 0.5 



(b) ($(-1/2), $(-2)) at r = 0.5 



Figure 2: {Ap/^ipix)) vs. ^ e (0 . . . r). The dots are TCSA results with cut 10 (48 states), the green 
dot-dashed and the solid black line is our result (4.31) up to the 3rd and 4th orders, respectively. The dotted 
red line is (4.33) and the violet dot-dot-dashed line is (4.31) without the term J22- The horizontal axes are 
positioned at the conformal (</5)q. 

4.2.2 Comparing the form factor expansion with TCSA 

Let us recall our main result that we will compare with the TCSA data: 

{0{x))l^ = {0)+Y. A, + . . . . (4.31) 

In the following we will refer to terms Dij with i + j = n as nth order terms. 

The work [13] addressed the calculation of the vacuum expectation value in large volume when 
the operator is close to one of the boundaries. With this setup they could neglect the effect of the 
other boundary and they worked up to the third order, so their formula in our notations is 

(O) + i?io + Ao + Ao . (4.32) 

One can try to improve this result for finite strips by considering also the analogous terms coming 
from the other boundary: 

3 

(0) + J](Ao + A*) • (4.33) 

i=l 

In order to show that our formula means a great improvement we will compare this expression with 
our result (4.31). 

Figure 2(a) shows the TCSA data with cut = 10 (48 states) together with our results for 
^M^/^(/9(,^)) for the boundary conditions (<&(6 = —2), <^(6 = —2)). In this case §0=9/3 = and only 
the even orders are non-zero. In particular, in this special case formula (4.33) coincides with our 
second order result. In Figure 2(b) we plot the case of the boundaries ($(6 = — 1/2),^>(6 = —2)). 
The TCSA has already converged, since changing the cut from A^ = 9 (38 states) to A^ = 10 (48 
states) the maximum change in the v.e.v. is ~ 0.0005. The plots clearly show that our formula gives 
much better results than (4.33) throughout the strip. We also plotted our result without the term 
J22 (violet dot-dot-dashed lines) in order to show that this non-trivial term is indeed present in the 
expansion (see the discussion in section 3.4.1). 

Figure 3 shows the v.e.v. for other pairs of boundary conditions. For the second pair changing 
the TCSA cut from A^ = 9 to A^ = 10 resuls in a change of order 10"^ in the v.e.v. For the (1, <I*(0)) 
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(a) (1,$(0)) at r = 3 



(b) ($(-1/2), $(-1/2)) at r = 2 



Figure 3: {Ap/^ipix)) vs. ^ S (0 . . .r). The dots are TCSA results with cut iV = 12 (45 states) in (a) and 
= 10 (48 states) in (b), the blue dashed, the green dot-dashed and the solid black line is our result (4.31) 
up to the 2nd, 3rd and 4th orders, respectively. The dotted red line is (4.33) and the violet dot-dot-dashed 
line is (4.31) with the sign of the term Dn flipped. The horizontal axes are positioned at the conformal ((^)q. 

boundaries raising the cut from = 8 (18 states) to = 12 (45 states) the change in the v.e.v. is 
only 0.05. 

In order to show that the conventions for the complex conjugation of the g factors are correct we 
plot also the result with the sign of the term Dn reversed (violet dashed lines). In Figure 3(b) we 
show also the second order result to show the convergence of the form factor expansion. It is clearly 
seen in both plots that adding the fourth order terms hardly changes the result, thus in these cases 
the form factor series is almost saturated by its first four orders. 

In Figure 4 we plot the one-point function for the (1, 1) boundaries for different strip widths. 
The TCSA is very convergent also in this case: a change in the cut from A^ = 8 (12 states) to 
A^ = 12 (29 states) causes a shift of order 10^^ in the v.e.v. The convergence of the form factor 
expansion is again obvious, however, for r = 2 our formula starts to deviate from the TCSA values, 
as a clear sign of the need for higher order terms. Our result is closer to the real values than the 
formula (4.33), but what is surprising is that the result of [13], expression (4.32), performs very well 
in the vicinity of the boundary. This can be understood by looking at the leading behaviour of the 
expectation value near the different boundaries [13]: 



(m2/V(x)) = (20^/^ '^'K (<^)*(0) + m''' '^K (1)^(0) + Oie'h , (4.34a) 



We have to check what happens when x is fixed and the volume R is changed. We know that the 
corrections to (4.32) are of the form exp(— Mi?), thus rescaling R is equivalent to a rescaling of M. 
Equation (4.34b) remains unchanged under such a rescaling in the leading order (and the next term 
is two orders higher in ^). In other words, the x-dependence of the v.e.v. agrees with the behaviour 
under a conformal rescaling, which means that for the l boundary the expectation value will not 
change much when the width of the strip changes. Obviously this is not true for (4.34a) 

Since our truncation of the form factor series is truly consistent only at x = R/2, in Figure 5 
we plot the deviation of the v.e.v. from the TCSA values at the middle of the strip for different 
volumes R. It is clearly seen on the logarithmic plot that for large enough volume (r > 2) the errors 
decrease exponentially with the volume and the exponents for the error of the consecutive orders 

^^We would like to thank Gerard Watts for discussing this question. 




(4.34b) 
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Figure 4: {hp/^tpix)) vs ^ G (0...r). The dots are TCSA results with cut iV = 12 (29 states), the blue 
dashed, green dot-dashed and the solid black line is our result (4.31) up to the 2nd, 3rd and 4th orders, 
respectively. The dotted red line is (4.33) and the asymmetric violet dot-dot-dashed line in (b) is (4.32). 
The horizontal axes are positioned at the conformal (</5)q. 

are approximately even-spaced. In particular, the sum of all four orders has the smallest and most 
rapidly decreasing error. A similar behaviour can be observed for other boundary conditions. 

5 Conclusions 

In this paper we considered vacuum expectation values of local operators in the presence of two 
integrable boundaries. We developed a consistent finite volume regularisation scheme to handle the 
various types of disconnected terms. The method presented in this work can be considered as a 
generalisation of the approach of [11] which dealt with expectation values with periodic boundary 
conditions. The key ingredients of our approach are the proper normalisation of the boundary state 
in a finite volume and the knowledge of the finite volume form factors. 

Apart from the task of evaluating disconnected terms of the form factors the two-boundary 
setting poses additional difficulties for gadp 7^ 0. In this case the boundary states include zero- 
momentum particles corresponding to the poles of the amplitudes Kj{6) at = 0. We showed 
that our finite volume evaluation scheme automatically provides a regularisation of the resulting 
singularities. We developed a new method to extract the finite parts (see Theorem 1 in section 2.2) 
which can be applied also in other cases, e.g. two-point functions at finite temperature. 

As a by-product of our formalism we were able to extract the second order terms of the ground 
state energy on the strip (subsection 2.2). This way we confirmed the boundary Thermodynamic 
Bethe Ansatz (BTBA) equations and we showed that they yield the correct second order terms 
even in the presence of zero-momentum particles [gaQp 7^ 0), for which no rigorous derivation of 
the BTBA is known. Also, we gave further independent support for the normalisation condition 
9i ~ which was proven by other means in [19, 20]. 

It is possible to compute higher order terms of the vacuum expectation value (0(x))^^ using the 
methods laid out in this work; in particular the finite volume normalisation of higher multi-particle 
contributions to the boundary state follows in a straightforward manner. However, the evaluation 
of the L — )• CO limit poses technical difficulties already at the next orders [Dij with i+j > 4). There 
will appear new types of singularities which are very similar to those encountered in the evaluation 
of thermal correlation functions (see subsection 3.4.3). These difficulties call for new summation 
schemes to separate the divergent and finite parts of the relevant finite volume series. Alternatively, 
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Figure 5: {M'^I^^{x = |)) - (m2/V(|)>tcsa ^-^^ ^) ^^^^ ^CSA cut TV = 12 (29 states) 
as a function of r. The violet dot-dot-dashed, the blue dashed, the green dot-dashed and the black 
solid line is our result (4.31) up to the 1st, 2nd, 3rd and 4th order, respectively. The red dotted line 
is formula (4.33). For r > 2 the errors decrease exponentially with r. 

it would be interesting to develop a general infinite volume regularisation scheme along the lines of 
[17, 6], which could be compared to the results of the present approach. 

There is plenty of room for future extensions of this work. First of all, one can consider expecta- 
tion values with respect to a boundary excited state. In this case one has to develop an appropriate 
modification of the finite volume boundary state formalism. A further interesting problem is the 
evaluation of (space and time dependent) correlation functions on the strip, which can be performed 
along the lines of the present paper. Future work is needed to generalise our methods to non-diagonal 
scattering theories which are relevant to a number of condensed matter problems [56]. Also, it needs 
to be clarified whether some of the present results also apply to massless scattering theories. In the 
absence of a mass gap the convergence of the spectral series is not guaranteed anymore and one 
does not have control over residual finite size corrections to form factors. However, in some cases it 
is possible to apply a modified version of the form factor program to massless theories [57, 58]. The 
study of a finite volume regularisation scheme for massless theories is left as an open problem. 
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Appendix 

A Regularisation of the divergent sums 

For a given L let us denote by 6i the solutions of the quantisation condition 



Q{e) = mLs\nh.e + (5(26*) = 27r/ 
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where 6{6) is the elastic phase shift defined by 

S{e) = -e'^'^^\ 5{-9) = -5{e). 

We also introduce the density of states 

(d) = '^^jp- = ml cosh 6 + 2ip{2e) , 
da 

where ip{e) = 5' (6). 

Theorem 1: Let f(9) be a symmetric function which apart from a double pole at 9 = is 
analytic in a neighbourhood of the real axis: 

f{9)^^ as 9^0. 



Then the expression 

G 



mL 



has a regular behaviour at large L with the L — )■ oo limit given by 



lim S{L) = If + Kf, 

L—^oo 



where 



Proof: 

We express the sum in the form of a complex integral for a finite L and perform the L — )■ oo limit 
afterwards. First of all, the summation over the rapidities can be replaced by a sum over contour 
integrals around the solutions of e**^ + 1 = 0: 



1 V ^ ^ 



2 ^ / 27r + 1 ' 



(A.2) 



where now the summation is over I G Z + 1/2. The integration contours can be transformed into 
two distinct curves: 

• the first starting from 9 = oo + ie running to e + ie, crossing the real axis between and the 
first solution of (A.l), then running from 9 = e — ie to 9 = oo — ie; 

• and a similar curve around the negative real axis, with the same counter-clockwise orientation. 

These two curves can be joined to form a single contour encircling the whole real axis. In doing this, 
one picks up the residue at ^ = 0, which is given by 

1 f(9) G G 

- Res .J^,:, ' — = -i—pi(0) = -i— (mL + 2(p(0)) . 

Therefore the 0{L) terms cancel and we obtain 
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Note that in the two integrals the L dependence is only contained in Q{9). We may now perform 
the L — )• oo limit for a fixed e > 0. To do this first of all note that the integrand has become 
bounded (we stay away from the real axis) therefore one may exchange the limit L — t- oo with the 
integration. Also, observe that 

lim .^.^ , . . = 1 , lim . , = . 

L^oo e^Qi^+^s) + 1 L^oo e«Q{^'-«e) + 1 

Putting everything together one obtains 

/co+ie 1 /~i 

The integral above can be pulled back to the real axis after an appropriate regularisation. We use 
the identity 

roo+ie _cosh0 



0= / —G rr- (A.4) 

J-oo+is smh 6 

which can be proven by deforming the contour to Im^ = m/2. Moreover, the function above has 
exactly the same singularity structure around = as the integrand in (A. 3). One may therefore 
subtract (A.4) from (A. 3) and pull back the integration contour to the real axis. The final result is 
thus given by 

B Properties of the four-particle form factor 

In this appendix we consider the behaviour of the function 

{61 + ^,-01+^,-62,62) (B.l) 

near 61,62 = 0. It follows from the form factor axioms that for infinitesimal 61^2 the function above 
is antisymmetric in both variables and disappears whenever one of the rapidities goes to zero whith 
the other one being kept fixed. Therefore the only allowed singularity is of the form 



(71(72 

However, it is easy to see that this pole does not appear. Setting 

di = 6 + e, 62 = 6, (B.2) 

and taking e — ?• results in [11] 

lim F^{6 + e + iTT,-6-£ + in, -6, 6) = (Fg{6, -6) - 2^{26)Fg) , (B.3) 

which is well-defined for every 6. Therefore a pole of cannot appear and it follows that 
lim Fg{6i + in, -61 + in, -62,62) = lim {Fg{6, -6) - 2^{26)Fg) = . 

61.O2 — ^0 6 — ^0 

Moreover, the first term of its Taylor expansion is proportional to 6162, as required in (3.37) to 
cancel the double pole of Ka{—6)Kp[6). 
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C Summary of our results 

In this appendix we collected the results of section 3. The vacuum expectation value is given by 

hi 

We calculated the terms with i + j < 4. The relative magnitudes of the individual contributions 
depend on x. However, a consistent ordering can be achieved by putting x = R/2. This way 
Dij behaves in the large R limit as 0{e~^^^^'^"^^^'^). A pictorial representation of the individual 
contributions can be found in Figs. 6 and 7. 

Terms of order e~™^/^ 

^10 - e 

Terms of order e~™^ 



D20 = \ j ^K^{e)F?{-e,e)e 



Ol a /IN 2m cosh 9 x 



Of a n\„~2mcosh6 (R—x) 



Terms of order e 3mi?/2 



1)30 = i|^A^„(0)|F3O(-e,^,O)e-™(2coshm). 

4 J 27r \ sinh 6 J 

ga f d0 f T?0(- a p\T^ fn\-2mcoshe(R-x)-mx 2(^/3 COsh 6 rn(2R~x)\ 

Du = - J - ^F,{^7^,-e,e)K,{e)e e >j 



37 



-<8) 



(a) Di 



(b) D20 



(d) L»4 



-<8>- 



(e) D^ 



(f) D31 



(g) I22 



Figure 6: Diagrams representing the different contributions to the v.e.v. Sohd hues correspond 
to pairs of particles with opposite momenta. Dashed hues correspond to the propagation of zero- 
momentum particles. The diagrams depicting the terms Dij with j > i can be obtained by switching 
the roles of the two boundaries. 
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(a) J22 



+ 



(b) D12 



+ 



(c) D2 



Figure 7: Diagrams representing contributions with singular pieces. Solid lines correspond to pairs of 
particles with opposite momentum. Dashed lines correspond to the propagation of zero-momentum 
particles. 
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Terms of order e 




D22 



I22 + J22 , where 




2m cosh 6ix—2m cosh 02 (R—x) 



and 



J22 = FO{in,0){ [^(KU-e)Ks{e)e-'"'^°^'^'^-^-^^^ 

[J ZTT \ 4 sinh 61/8 J 



References 

[1] P. Calabrese and J. L. Cardy, "Time-dependence of correlation functions following a quantum 
quench," Phys. Rev. Lett. 96 (2006) 136801, arXiv: cond-niat/0601225. 

[2] V. Gritsev, E. Demler, M. Lukin, and A. Polkovnikov, "Analysis of quench dynamics of 

coupled one dimensional condensates using quantum sine gordon model," Phys. Rev. Lett. 99 
(2007) 200404, arxiv: cond-niat/0702343. 

[3] D. Fioretto and G. Mussardo, "Quantum quenches in integrable field theories," 
arxiv: 0911 . 3345 [cond-mat . stat-mech] . 

[4] A. LeClair, F. Lesage, S. Sachdev, and H. Saleur, "Finite temperature correlations in the 
one-dimensional quantum ising model," Nuclear Physics B 482 (1996) 579, 
arxiv: cond-mat/9606104. 

[5] R. Konik, A. LeClair, and G. Mussardo, "On Ising correlation functions with boundary 
magnetic field," Int. J. Mod. Phys. All (1996) 2765-2782, arXiv:hep-th/9508099. 

[6] D. Schuricht and F. H. L. Essler, "Dynamical response functions in the quantum ising chain 
with a boundary," J. Stat. Mech.: Theor. Exp. (2007) P11004, 
arXiv : 0709 . 1809 [cond-mat . str-el] . 

[7] A. LeClair and G. Mussardo, "Finite temperature correlation functions in integrable QFT," 
Nucl. Phys. B552 (1999) 624-642, arXiv:hep-th/9902075. 

[8] O. A. Castro-Alvaredo and A. Fring, "Finite temperature correlation functions from form 
factors," Nucl. Phys. B636 (2002) 611-631, arXiv:hep-th/0203130. 

[9] H. Saleur, "A comment on finite temperature correlations in integrable QFT," 
Nucl. Phys. B567 (2000) 602-610, arXiv:hep-th/9909019. 



40 



[10] G. Mussardo, "On the finite temperature formalism in integrable quantum field theories," 
J. Phys. A34 (2001) 7399-7410, arXiv:hep-th/0103214. 

[11] B. Pozsgay and G. Takacs, "Form factors in finite volume ILdisconnected terms and finite 
temperature correlators," Nucl. Phys. B788 (2008) 209-251, arXiv : 0706 . 3605 [hep-th] . 

[12] G. Takacs, "Finite temperature expectation values of boundary operators," 
Nucl. Phys. B805 (2008) 391-417, arXiv: 0804.4096 [hep-th]. 

[13] P. Dorey, M. Pillin, R. Tateo, and G. M. T. Watts, "One-point functions in perturbed 
boundary conformal field theories," Nucl. Phys. B594 (2001) 625-659, 
arXiv : hep-th/0007077. 

[14] L. Samaj and Z. Bajnok, "Exactly solvable model of the 2d electrical double layer," 
Physical Review E 72 (2005) 061503, arXiv: coiid-mat/0508410. 

[15] B. Pozsgay and G. Takacs, "Form factors in finite volume I: form factor bootstrap and 

truncated conformal space," Nucl. Phys. B788 (2008) 167-208, arXiv: 0706 . 1445 [hep-th]. 

[16] F. H. L. Essler and R. M. Konik, "Finite-temperature lineshapes in gapped quantum spin 
chains," Phys. Rev. B78 (2008) 100403, arXiv: 0711 . 2524 [cond-mat . str-el] . 

[17] F. H. L. Essler and R. M. Konik, "Finite temperature dynamical correlations in massive 
integrable quantum field theories," J. Stat. Mech. (2009) P09018, 
arXiv: 0907. 0779 [cond-mat]. 

[18] S. Ghoshal and A. B. Zamolodchikov, "Boundary S matrix and boundary state in 

two-dimensional integrable quantum field theory," Int. J. Mod. Phys. A9 (1994) 3841-3886, 
arXiv : hep-th/9306002. 

[19] Z. Bajnok, L. Palla, and G. Takacs, "Finite size effects in quantum field theories with 

boundary from scattering data," Nucl. Phys. B716 (2005) 519-542, arXiv:hep-th/0412192. 

[20] Z. Bajnok, L. Palla, and G. Takacs, "Boundary one-point function, Casimir energy and 
boundary state formahsm in D+1 dimensional QFT," Nucl. Phys. B772 (2007) 290-322, 
arXiv : hep-th/0611 176. 

[21] F. A. Smirnov, "Form-factors in completely integrable models of quantum field theory," Adv. 
Ser. Math. Phys. 14 (1992) 1-208. 

[22] B. Berg, M. Karowski, and P. Weisz, "Construction of Green Functions from an Exact S 
Matrix," Phys. Rev. D19 (1979) 2477. 

[23] A. B. Zamolodchikov, "Two point correlation function in scaling Lee- Yang model," 
Nucl. Phys. B348 (1991) 619-641. 

[24] G. Delfino and G. Mussardo, "The Spin spin correlation function in the two-dimensional Ising 
model in a magnetic field at T = T(c)," Nucl. Phys. B455 (1995) 724-758, 
arXiv : hep-th/9507010. 

[25] G. Delfino and P. Simonetti, "Correlation Functions in the Two-dimensional Ising Model in a 
Magnetic Field at T = T^," Phys. Lett. B383 (1996) 450-456, arXiv:hep-th/9605065. 



41 



[26] I. Affleck and A. W. W. Ludwig, "Universal noninteger 'ground state degeneracy' in critical 
quantum systems," Phys. Rev. Lett. 67 (1991) 161-164. 

[27] P. Dorey, D. Fioravanti, C. Rim, and R. Tateo, "Integrable quantum field theory with 
boundaries: The exact g-function," Nucl. Phys. B696 (2004) 445-467, 
arXiv : hep-th/0404014. 

[28] P. Dorey, C. Rim, and R. Tateo, "Exact g-function flow between conformal field theories," 
arXiv: 0911. 4969 [hep-th] . 

[29] M. Liischer, "Volume Dependence of the Energy Spectrum in Massive Quantum Field 
Theories. 1. Stable Particle States," Commun. Math. Phys. 104 (1986) 177. 

[30] T. R. Klassen and E. Melzer, "On the relation between scattering amplitudes and finite size 
mass corrections in QFT," Nucl. Phys. B362 (1991) 329-388. 

[31] Z. Bajnok and R. A. Janik, "Four-loop perturbative Konishi from strings and finite size effects 
for multiparticle states," Nucl. Phys. B807 (2009) 625-650, arXiv : 0807 . 0399 [hep-th]. 

[32] Y. Hatsuda and R. Suzuki, "Finite-Size Effects for Multi-Magnon States," 
JHEP 09 (2008) 025, arXiv : 0807 . 0643 [hep-th]. 

[33] B. Pozsgay, "Luscher's mu-term and finite volume bootstrap principle for scattering states 
and form factors," Nucl. Phys. B802 (2008) 435-457, arXiv: 0803 .4445 [hep-th]. 

[34] A. LeClair, G. Mussardo, H. Saleur, and S. Skorik, "Boundary energy and boundary states in 
integrable quantum field theories," Nucl. Phys. B453 (1995) 581-618, 
arXiv : hep-th/9503227. 

[35] P. Dorey, A. Pocklington, R. Tateo, and G. Watts, "TBA and TCSA with boundaries and 
excited states," Nucl. Phys. B525 (1998) 641-663, arXiv:hep-th/9712197. 

[36] Z. Bajnok, C. Rim, and A. Zamolodchikov, "Sinh-Gordon Boundary TBA and Boundary 
Liouville Reflection Amplitude," Nucl. Phys. B796 (2008) 622-650, 
arXiv: 0710. 4789 [hep-th]. 

[37] G. Takacs, "Form factor perturbation theory from finite volume," 
Nucl. Phys. B825 (2010) 466-481, arXiv: 0907. 2109 [hep-th]. 

[38] B. Pozsgay, "work in progress,". 

[39] J. L. Cardy, "Conformal invariance and the Yang-Lee edge singularity in two-dimensions," 
Phys. Rev. Lett. 54 (1985) 1354-1356. 

[40] J. L. Cardy, "Boundary Conditions, Fusion Rules and the Verlinde Formula," 
Nucl. Phys. B324 (1989) 581. 

[41] J. L. Cardy and D. C. Lewellen, "Bulk and boundary operators in conformal field theory," 
Phys. Lett. B259 (1991) 274-278. 

[42] D. C. Lewellen, "Sewing constraints for conformal field theories on surfaces with boundaries," 
Nucl. Phys. B372 (1992) 654-682. 

[43] I. Runkel, "Boundary structure constants for the A-series Virasoro minimal models," 
Nucl. Phys. B549 (1999) 563-578, arXiv:hep-th/9811178. 



42 



[44] J. L. Cardy and G. Mussardo, "S Matrix of the Yang-Lee Edge Singularity in Two- 
Dimensions," Phys. Lett B225 (1989) 275. 

[45] A. B. Zamolodchikov, "Thermodynamic Bethe Ansatz in relativistic models. Scaling three 
state Potts and Lee- Yang models," Nucl. Phys. B342 (1990) 695-720. 

[46] A. B. Zamolodchikov, "Mass scale in the sine-Gordon model and its reductions," 
Int. J. Mod. Phys. AlO (1995) 1125-1150. 

[47] F. A. Smirnov, "The perturbed c<l conformal field theories as reductions of the Sine-Gordon 
model," Int. J. Mod. Phys. A4 (1989) 4213-4220. 

[48] F. A. Smirnov, "Reductions of the sine-Gordon model as a perturbation of minimal models of 
conformal field theory," Nucl. Phys. B337 (1990) 156-180. 

[49] A. B. Zamolodchikov, "Two point correlation function in scaling Lee- Yang model," 
Nucl. Phys. B348 (1991) 619-641. 

[50] V. Fateev, S. L. Lukyanov, A. B. Zamolodchikov, and A. B. Zamolodchikov, "Expectation 
values of local fields in Bullough-Dodd model and integrable perturbed conformal field 
theories," Nucl. Phys. B516 (1998) 652-674, arXiv:hep-th/9709034. 

[51] P. Dorey, I. Runkel, R. Tateo, and G. Watts, "g-function flow in perturbed boundary 
conformal field theories," Nucl. Phys. B578 (2000) 85-122, arXiv:hep-th/9909216. 

[52] V. P. Yurov and A. B. Zamolodchikov, "Truncated Conformal Space Approach to Scaling 
Lee- Yang model," Int. J. Mod. Phys. A5 (1990) 3221-3246. 

[53] G. Feverati, K. Graham, P. A. Pearce, G. Z. Toth, and G. Watts, "A renormalisation group 
for TCSA," arXiv:hep-th/0612203. 

[54] M. Kormos, I. Runkel, and G. M. T. Watts, "Defect flows in minimal models," 
JHEP 11 (2009) 057, arXiv : 0907 . 1497 [hep-th] . 

[55] I. Runkel, "Phd thesis." http://www.nith.kcl.ac.uk/staff/i_runkel/PDF/phd.pdf. 

[56] F. H. L. Essler and R. M. Konik, "Applications of massive integrable quantum field theories 
to problems in condensed matter physics." http://arxiv.org/abs/cond-niat/0412421, 2004. 

[57] G. Delfino, G. Mussardo, and P. Simonetti, "Correlation functions along a massless flow," 
Phys. Rev. D51 (1995) 6620-6624, arXiv:hep-th/9410117. 

[58] D. C. Cabra, A. Honecker, G. Mussardo, and P. Pujol, "A non-perturbative approach to the 
random-bond Ising model," J. Phys. A30 (1997) 8415-8426, arXiv :hep-th/9705180. 



43 



